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^ Abstract 

^ We investigate the Zel'dovich effect in the context of ultra-cold, harmonically trapped quantum 



(N 



gases. We suggest that currently available experimental techniques in cold-atoms research offer 
an exciting opportunity for a direct observation of the Zel'dovich effect without the difficulties 



^ imposed by conventional condensed matter and nuclear physics studies. We also demonstrate an 

I interesting scaling symmetry in the level rearragements which has heretofore gone unnoticed. 
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I. INTRODUCTION 

The Zel'dovich effect (ZE)P occurs in any quantum two-body system for which the con- 
stituent particles are under the influence of a long range attractive potential, suplemented by 
a short-range attractive two-body interaction, which dominates at short distances. The sys- 
tem flrst found to exhibit the ZE consists of an electron experiencing an attractive long-range 
Coulomb potential, which at short distances, is modifled by a short-range interaction.'^ 

The characterizing feature of the ZE in this scenario is that as the strength of the at- 
tractive two-body interaction reaches a critical value {i.e., when a two-body bound state is 
supported in the short-range potential alone) , the S^-wave spectrum of the distorted Coulomb 
problem evolves such that the ground state IS level plunges down to large negative energies 
while simultaneously, the flrst radially excited 25^ state rapidly falls to flU in the "hole" left 
by the ground state level. This also occurs for higher levels in which, generally, the (n + 1)5' 
level replaces the nS level. This so-called level rearrangement is the signature of the ZE, 
and continues as the strength of the two-body interaction is further increased to support 
additional low-energy scattering resonances (see e.g.. Fig. 2 of Ref. jT]). 

Recently, Combescure et. al have revisited the ZE in the context of "exotic atoms" ^^ 
where a negatively charged hadron replaces the electron, and the short-range interaction is 
provided by the strong nuclear force. However, tuning the short-range interaction in exotic 
atoms implies that one must be able to adjust the nuclear force in the laboratory, which 
is a formidable task. Indeed, while it is theoretically easy to adjust the strength of the 
short-range interaction between the particles in any of the systems above, the experimental 
reality is very diflFerent. As a result, direct experimental observation of the ZE has been 
lacking, in spite of suggestions for its observation in quantum dots,'^' Rydberg atoms,Hl and 
atoms in strong magnetic flelds.'^ 

In this article, we explore the possibility for a direct observation of the ZE in harmonically 
trapped, charge neutral, ultra-cold atomic gases. The charge neutrality of the atoms ensures 
that the long-range attractive potential is provided solely by the isotropic harmonic oscillator 
trap, while the short-range two-body interaction is naturally present owing to the two-body 
5-wave scattering, which is known to dominate at ultra-cold temperatures. Moreover, the 
short-range interaction between the atoms is completely tuneable in the laboratory via the 
Feshbach resonance.^ The multi-channel Feshbach resonance can be treated in a simpler 



single-channel model by a finite-range, attractive two-body interaction, supporting scatter- 
ing resonances. Thus ultra-cold atoms, at least in principle, provide all of the necessary 
ingredients for the experimental observation of the Zel'dovich effect. 

The plan for the remainder of this paper is as follows. In Section II, we establish a 
deep connection between the level rearrangements and the two-body energy spectrum as 
characterized by the 5-wave scattering length, a. This connection allows us to make contact 
with recent experimental results on ultra-cold two-body systems,^ from which we suggest 
that a direct observation of the ZE is possible. Then, in Section III, we investigate the 
influence of the range of the two-body interaction on the level rearrangements examined 
in Section II. In particular, we reveal an interesting scaling symmetry of the two-body 
energy spectrum which has not been noticed before. Finally, in Section IV, we present our 
concluding remarks. 

II. ZEL'DOVICH EFFECT IN ULTRA-COLD ATOMS 

A. Universal two-body energy spectrum 

The two-body spectrum for a pair of harmonically trapped ultra-cold atoms is obtained 
from the following Hamiltonian, 

^ = ll + ^ + ^^"^''^ + ^^"^''^ + ^^^^'"^ " '^]^ ' ^^^ 

where each atom has a mass of M and Ksi^d^i — r2|) is a short-range potential. Introducing 
the usual relative, r = ri — r2, and centre of mass, R = (ri + r2)/2 coordinates, and noting 
that the centre of mass motion may be separated out, the associated Schrodinger equation 
in the 5- wave channel reads^ 



- ^u"{r) + ^MujWir) + VsR{r)u{r) + ^i^^^Jluir) = Eu{r) , (2) 

where u{r) = r^^-^y^ip^r) is the reduced radial two-body wave function, primes denote 
derivatives, d is the dimension of the space, and E is the relative energy of the two-body 
system. Defining the dimensionless variables 77 = 2E/huj^ lose = \/h/Muj and x = r/\/2£osc^ 
Eq. (2) may be written as 

— u^\x) + x'^u(x) + Vsr(x)u(x) -\ u(x) — Tjuix) = 0, (3) 



where Vsr{x) = 2Vsr/Jtu. 

Exact analytical solutions to (JSJ) exist \/d if the potential is taken to be an appropriately 
regularized zero-range contact interaction. For d = 3 the spectrum is described by^'^ 

a T{l/4-E/{2nuj)) 



for d = 1 we haveP 



Lsc ^/2T{3/4-E/{2hoj))' 
^r(3/4 - E/{2hu)) 



(4) 



^OSC 



and for d = 2^ 



a r{l/i-E/{2huj)) ' ^^^ 

V^(l/2 - E/(2hw)) = In -^ + 2 In 2 - 27 . (6) 

2a^ 

In the above, a is the 5-wave scattering length in Vsr alone, r(-) is the gamma function, 
ip{-) is the digamma function and 7 = 0.577215665... is the Euler constant. ^^ Note that in 
any dimension, the two-body spectrum is universal in the sense that the relative energy, £", 
is determined entirely by the scattering length. Thus, even for a two-body potential with 
finite range, 6, it has been shown that provided b <C iosc (practically speaking, b/£osc ^ 0.01), 
the same two-body energy spectrum as described above will be obtained for an arbitrary 
two-body interaction evaluated at the same scattering length.'^ 

B. Level rearrangements 

In this section, the level rearrangements {i.e., the ZE) exhibited by the two-body energy 
spectrum are investigated. The results presented here are strictly for three-dimensions (3D), 
although analogous findings are also observed in other dimensions. We will focus on three 
diflFerent interaction potentials, viz., a finite square well (FSW), the modified Poshl- Teller 
potential^^ and an exponential potential^ 



VsR{r) = { 



-VoOib-r) 

-Vosech\r/b) (7) 

-Voexp(-r/6) , 



respectively. In the above, G(-) is the Heaviside step function, and Vq is the depth of the 
potential. The 3D s-wave scattering lengths for the potentials are given by (in the same 
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FIG. 1: Left panel: 5- wave two-body energy spectrum versus strength for all 3 model potentials 
with fixed b/£osc — 0.01. Open circles, triangles and squares represent the numerical integration 
of Eq. (3) for the FSW, Poshl- Teller and exponential potentials, respectively. The vertical dashed 
line represents the critical strength, gc, of all three potentials while the horizontal dashed lines 
are energy values at a = zboc. The strength axis is scaled so that the critical strength values lie 
along the same vertical dashed line. Right panel: Energy versus scattering length for all 3 model 
potentials. The same symbols as the plots in the left panel are used. The vertical dot-dash line 
indicates a = 0. In all 4 plots, the solid black line is the exact expression obtained from Eq. Q. 
Units are scaled as discussed in the text. 



order as the potentials listed above )'^2'I2HIS1 
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6 1 



tan^ 
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6 7 + V^(A) + |cot^A/2 



(8) 






b h-f + lng 

where g = MVob'^/h'^ is the dimensionless strength of the potential, A = (1 — vT^F4^)/2, 
and Jo(') and lo(') are the zeroth order Bessel functions of the first and second kind, re- 
spectively.'^ 



We proceed by numerically integrating Eq. ^ for each of the three potentials. For our 
numerics, we have set h = uj = 1 and M = 2 to be consistent with the numerical results 
in Refs. pfTOj . We plot our numerical results (open symbols) for the relative energy, E (in 
units of hu))^ as characterized by both the strength, g^ and the 5- wave scattering length, a, in 
Figure 1. The left panel illustrates the level rearrangements as the strength, g^ is increased 
beyond the first scattering resonance, whereas the right panel illustrates the relative energy, 
£", as determined by the 5-wave scattering length. The level rearrangements shown in the 
left panels illustrate the 2S level replacing the IS level at the first scattering resonance, 
while the IS level dives down to large negative values. 

A further examination of Fig. 1 reveals that while h/losc = 0.01 for all three potentials, the 
level rearrangements displayed in the left panels exhibit noticeable diflFerences. In particular, 
we see that the FSW has a much sharper drop ai g = gc^ than the Poshl- Teller or exponential 
potentials. These level repulsions, or "anticrossings" , are known to be as a result of the levels 
belonging to the same SO (2) symmetry of the Hamiltonian, while the mixing of the levels 
is dependent on how rapidly the short-range potential "shuts off" P^ 

The underlying message here is as follows. While all three plots on the left of Fig. 1 
display the Zel'dovich effect, namely they all undergo level rearrangement at some value of 
the strength parameter g^ all three different potentials map on to the same E vs. a curve, 
as illustrated in the right panel of Figure 1. This reaffirms that while the details of the 
ZE are sensitive to the form of the two-body interaction, the energy dependence on the 
scattering length, a, is indeed universal. It is also worthwhile pointing out that the solid 
curves in the left panels of Fig. 1 are obtained from substituting the expressions for the 
scattering length, Eq. (8), into the Eq. (4), which is exact only for a zero-range interaction. 
However, it is clear that the numerically obtained open symbols closely follow the solid curve 
derived from Equation (4). Thus, for b/£osc <^ 1, the level rearrangements in harmonically 
trapped two-body systems interacting via a finite, short-range potential, are all equivalent 
to a zero-range interaction. Viewed another way, given a set of data for E vs. a, there 
must exist some quantum two-body system (z.e.^ the two-body potential need not be known 
explicitly) whose E vs. g dependence exhibits the Zel'dovich effect. This observation has 
some interesting implications, which we further explore in the next subsection. 



C. Flow of the Spectrum 



In order to make the connection between the E vs. g and E vs. a curves more apparent, 
we now study the "flow" of the two-body energy spectrum. Although we focus our attention 
to the FSW, the same analysis holds for any other potential. 
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FIG. 2: (Color online) Flow of the spectrum for the FSW. Boxed and double arrows (green online) 
represent the flow of the flrst excited state while circles and single arrows (red online) follow the flow 
of the ground state. Identical points in each panel are labeled by the same letter. For example, the 
point B on the left is the exact same data point as B on the right but subject to the transformation 
in Equation ([s]). The horizontal dotted lines in both panels correspond to the asymptotic values 
for the energy E at \a\ = oo, as in Figure 1. See also the discussion preceding Eq. (19) in Section 
III. 



In the left panel of Fig. 2, we note that as g is increased from zero, the energy only slightly 
varies from the unperturbed energy, until the critical strength, gc^ is reached at which point 
the Zel'dovich effect occurs. In the right panel of Fig. 2, the same flow is illustrated, but 
this time in terms of the scattering length. The lower flow in the left panel (red online) 
illustrates that the trajectory of the ground state A^B^C^Dis continuous though 
the resonance Sit g = gc- However, as we follow the same path in £" vs. a, the point B flows 
out to a ^ — oc while C and D flow in from a ^ +oo and then to a ^ 0. Thus, while 
the flow for the energy spectrum in ^'-space is continuous, the flow in a-space appears to 
be disconnected. Similarly for the flrst excited state (green online) where the B^ and C^ 
flow is continuous in g'-space, but rapidly branches oS to a ^ —oo and a ^ 0, in a-space. 



respectively. The continuous flow in ^'-space suggests that the a-space spectrum is more 
appropriately viewed on the topology of a cylinder, where a = ±oc may be identified. 

1. Cylindrical Mapping 

The observations made above suggest that we map the E vs. a spectrum onto the surface 
of a cylinder. The details of this mapping are closely related to the mapping of the real line 
(in our case, the scattering length) onto the unit circle, S'^, followed by constructing the 
Cartesian product, M x S''^, with M identified with the energy, E. The essential point of this 
mapping is to provide a more natural interpretation for the two-body E vs. a spectrum. 

To this end. Fig. 3 illustrates a series of "snapshots" which show how the the original 
E vs. a spectrum is mapped onto the surface of a cylinder. Each of the 6 panels in Fig. 3 



^^^ 


J^il'l^ 




/^ 












FIG. 3: The E vs. a spectrum being rolled onto a cylinder. Top from left to right: a = —10 to 
a = 10, a = —30 to a = 30, a = —35 to a = 35. Bottom from left to right: a = —40 to a = 40, 
a = —90 to a = 90, a = — oc to a = oo. The thick vertical line (red online) represents a = 0. The 
symmtery axis of the cylinder is the E axis while the azimuthal angle is connected to the scattering 
length, a. 



should be viewed as an intermediate step in taking the E vs. a spectrum and roUing it onto 
a cyhnder. In Fig. 4, we present the complete mapping of the E vs. a spectrum up to the 
4-th excited state of the bare harmonic trap. This "Zel'dovich spiral" (ZS) may now be 
explicitly connected to the level rearrangements discussed in the left panel of Fig. 1 above. 
Indeed, we observe that the flow of the spectrum shown in the left panels of Fig. 2 
correspond to clockwise (CW) rotations about the Zel'dovich spiral. That is, increasing the 
strength of the two-body interaction corresponds to moving along the ZS in a CW direction, 
with the starting point {i.e., the front of the cylinder) along the thick vertical line (red 
online) in Figure 4. A CW rotation of n puts us on the back of the cylinder, or a = — oc. 




FIG. 4: The complete mapping of the E vs. a spectrum onto the surface of a cylinder illustrating 
the "Zel'dovich spiral". The solid vertical line (red online) identifies the unperturbed system for 
which a = 0. Every 27r winding along the spiral corresponds to a complete level rearrangement 
e.g., 2S -^ IS after 2ti rotations. The lower solid circle indicates the unperturbed IS level, whereas 
the upper solid circle corresponds to the 2S level. 



whereas a counter-CW rotation of tt takes us to a = +oc {i.e., the azimuthal angle \(j)\ = n 



is a branch point). 

To see how the ZS naturaUy contains the level rearrangements, let us first begin at 
E = 3/2, which in Fig. 4 is represented by the lower solid circle along the vertical line. As 
we move in a CW rotation along the spiral, the E = 3/2 (a = 0, = 0) goes to £" ^ 1/2+ 
for large negative values of a, and finally to £" = 1/2 at a = — oc (lower dotted curve in 
the right panel of Figure 2). A further infinitesimal rotation takes us to £" = l/2~ at large 
positive values of a, and finally to £" ^ — oc at a = after a full 27r rotation; we have just 
followed the fiow of the IS level in the left panel of Fig. 1, viz., A^B^C^D^---. 
Similarly, the upper solid circle in Fig. 4 corresponds to £" = 7/2, which as we rotate CW, 
evolves to £" ^ 5/2+ for large negative values of a, £" = 5/2 at |a| = oc (upper dotted curve 
in the right panel of Figure 2), and subsequently to £" = 3/2 at = 27r; this description is 
precisely the fiow of the 2S level in the left panel of Figure 2. If we were to the continue 
with our CW rotation {i.e., continue increasing the strength, g)^ we would then evolve from 
E = 3/2^ 1/2 at = 37r followed hy E ^ -oc at = An. 

In our opinion, viewing level rearrangements in this way is more natural than the original 
E vs. a spectrum in M^. We see that critical strengths, gc^ correspond to CW rotations of 
odd multiples of n whereas a complete level rearrangement occurs for even multiples of n. 
In general, the (n + 1)5' level, with En = 2n + 3/2 (n = 0, 1, 2, ...) will eventually evolve to 
the IS level after 2n7T CW rotations along the spiral. 

The Zel'dovich spiral also helps to clarify several misconceptions about E vs. a spectrum 
in the literature. The spectrum is typically understood by taking a = ±oc separately, and 
assigning diflFerent interpretations to a ^ 0+ and a ^ 0~ . An example of this is a recent 
contribution is by Shea et a/,^^ where the authors describe the spectrum by first "starting 
from the far left" and making the interaction weaker and weaker as a ^ 0~ and then 
independently "starting from the right" and making the interaction stronger and stronger as 
a ^ 0+. On the ZS, nothing is ambiguous, since one always moves in a CW rotation along 
the spiral, corresponding to increasing the strength, g^ of the interaction; a complete level 
rearrangement occurs after we undergo an even muliple of n CW rotations. Furthermore, 
the "counter-intuitive" properties of the E V8. a spectrum discussed in Ref. [9j are now seen 
to be nothing more than a manifestation of the onset of the Zel'dovich effect. We find it 
rather surprising that the ZE has been present in the two-body E vs. a spectrum all along, 
but until now, has gone unnoticed. 

10 



2. Experimental Observations 

In a recent work, Stoferle et. alj have experimentally measured the binding energy as 
a function of the 5- wave scattering length between two interacting particles in a harmonic 
trap. This experiment highlights the versatility of trapped, ultra-cold atomic systems, in 
which an analytically solvable model, once only the purview of theoretical physics, has 
now been realized in the laboratory. Remarkably, the experimental results for the E vs. a 
spectrum are in excellent agreement with theory, (see Fig. 2 in Ref. [7J), even though the 
two-body interaction in the experiments is most certainly not a zero range interaction. 
Thus, the theoretical prediction that the E vs. a spectrum is universal has been confirmed 
experimentally. 

What has not been appreciated until now, however, is that the experimental E vs. a 
spectrum obtained in Ref. [7J is exactly equivalent to obtaining the ground state branch in 
the left panel of Fig. 2 (single arrows, red online). In other words, the work of Stoferle et. 
al^ has already been a direct experimental observation of the ground state branch of the two- 
body system exhibiting the Zel'dovich effect. We therefore suggest that further experiments 
along the lines of Ref. [7J be performed so that data corresponding to the double arrows and 
primed letters (green branch online) in the right panel of Fig. 2 may be obtained. If such an 
extension to the experiments in Ref. [7J is viable, then according to our analysis, this data 
would be exactly equivalent to the 2S branch (double arrows, green online, in the left panel 
of Fig. 2) undergoing the Zel'dovich eflFect. Therefore, just a few additional data points in 
the E vs. a spectrum, would provide for a direct experimental confirmation of the ZE for 
two interacting particles confined in a harmonic trap. 

III. LEVEL REARRANGEMENTS IN THE ZERO-RANGE LIMIT 

We close this work with a discussion of an interesting scaling symmetry present in the 
the level rearrangements. Specifically, we show that in the 6^0 limit, the entire two-body 
E V8. g spectrum is determined by only the first level rearrangement. 

Figure 5 illustrates the spectrum through the first three low-energy scattering resonances, 
viz.^ g = go^gi^g2 for the FSW. From this figure, we immediately notice the similarities 
between the level rearrangements as we move from one region to the next, along a fixed 
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FIG. 5: Several level rearrangements for a FSW interaction with h/losc — 0.01. Dotted horizontal 
lines (green online) are unperturbed energy values while dotted vertical lines (red online) are values 
at which a = 0. Each subsequent region corresponds to a CW rotation of 27r along the Zel'dovich 
spiral. The solid black circle is a representative data point which we wish to map back into the g^ 
region, as schematically illustrated by the open circle in ^o- Single arrows on the ^-axis indictate 
the various critical values, ^^n? in the n^^ region, while double arrows indicate the ^q for which 
the scattering length vanishes. 

value for the energy, E. Each region begins at a = 0, undergoes a level rearrangement, and 
then returns to a = 0. In the language of the ZS, each region corresponds to one complete 
27r CW rotation along the spiral. Figure 5 suggests that it may be possible to map every 
Qn {n 7^ 0) region onto g^) by some appropriate scaling of the ^'-axis. This mapping should 
ensure that a = in, say, gi^ matches with a = in 5^0 , and that the critical gi value in the 
first region overlaps with the critical g^ value in the zeroth-region, and so on. 



(n) 



Let us first define some useful nomenclature. We define g\^ as the n value of g for 
which a = (double arrows in Fig. 5). Next, gc^^ is defined as the n^^ critical g value; that 
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is the g value at which the n^^ level rearrangement occurs (single arrows in Fig. 5). Lastly, 
g is the value of g outside the region go which we intend to map back into go. For example, 
consider the point labeled by the solid dot in the g2 region of Figure 5. Considering this 
point, which we wish to map back into go (represented by the open circle in Fig. 5), we have 
g = 75, 5'c,2 = 257r^/4 and go — 59.67951594410.... The mapping that takes this g back 
into go is 

g = ^^~^Q )^^-Q ^ 18.84894603..., (9) 

9c,2 - go 

where gc^o = 7r^/4 is the zeroth critical g value. We may generalize this example to any 
region by employing the following prescription 



9c,n Qo 



{n^Q). (10) 



With this remapped value oi g^ we also have the associated energy, E. If the mapping is 
indeed exact, the energy E of the remapped point should be identical to the the energy for 
the same g value in go- 

In Fig. 6 we study this mapping for the FSW with h/losc = 0.01. At first glance, the 
left panel Fig. 6 appears to be show that the mapping is exact, but a closer examination 
of the spectrum for values of g near the resonance (right panel in Fig. 6) reveals noticeable 
discrepancies between data in the go and gn > go regions. Remarkably, even for b/£osc = 0.01, 
the mapping of the data from gi and g2 into go agree almost perfectly (z.e.^ the dashed (red 
online) and the dot-dashed (green online) curves, respectively). Regardless, the mapping 
given by Eq. (10) is not exact for any finite range, b. 

We can, however, show that this mapping becomes exact in the 6^0 limit by considering 
Equation Q. We write this expression in the notationally convenient form a{g) = R{E) 
where R{E) = Y{^/l-'E/{2nw\) ^^^ ^(^) ^ '^^{g)/^osc- Our goal is now to map g values in two 
different regions, gk and gk'^ onto a value in the region go and investigate the difference in 
their energy values. We denote 5E = E^' — E^ and note that we have the two expressions 
for the spectrum d{gk) = R{Ek) and d{gkf) = R{Ek'). The difference between these two 
expressions is d{gk') — d{gk) = R{Ek') — R{Ek) = R{SE + E^) — R{Ek). Taylor expanding 
up to first order in 5E gives 
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FIG. 6: Left: Equation (10) as applied to data from the FSW to different regions with h/losc — 0.01. 
Right: The same data as in the left-panel, but with g values surrounding the critical value, ^co- 
in both panels, the solid (blue online), dashed (red online) and dot-dashed (green online) curves 
correspond to the g^^gi and ^2 regions, respectively. In these figures, the remapped data from gi 
and g2 are indistinguishable on the scale of the plots. 



where 5a = a{gk') — o.{gk) and R'{Ek) = L . From Eq. (4) we may re-express Eq. (11) as 



5E = 



5a dEk 



^'{gk) dgk ' 

which upon noting that E^ = R~^{a{gk))^ becomes the implicit expression 

5a 



R'{a{gk)) = R 



5E 



(12) 



(13) 



Assuming SE to be small compared to Sa, we seek an asymptotic expression for 



R 



Sa \ V2T{l/4 - ^) 



: = 



(14) 



^SEJ r(3/4-^) • 

An application of Euler's reflection formulaP'and Stirling's approximation to the above gives 
the approximate expression 



^ , 5a\ 5E f Ti5a Stt 



(15) 



Equation (15), along with Eq. (13) gives 



5a 2 , I 5E 1 
5E- 7^ I V <5a R'{a{gk)) 



3 

2' 



(16) 
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With the approximation cot ^{x) c^ 7r/2, x <C 1 the difference in the energies becomes 

5E = ha, (17) 

5 

or, defining a{g) = bA{g), 

SE = ^ {A{g,>) - A{g,)) . (18) 



Equation (18) analyticaUy shows that SE ^ as 6 ^ 0, and the g values in the two diflFerent 
regions get mapped back into the the go region at the exact same energy. Therefore, the 
mapping of all subsequent regions, gn {n ^ 0) back onto g^ is exact in the zero-range 
limit. It is important to note that our analysis has not relied upon specifying the details 
of interaction, and so is equally valid for any short-range two-body interaction supporting 
bound states. 

It is also instructive to consider how the shape of the level rearrangements curves evolve 
as 6 ^ 0. The shape-dependence of the curves can be established by expanding the right 
hand side of Eq. ^ about g = gc (some resonant strength value) and the left hand side 

about E = Ec where, in 3D, Ec = 1/2,5/2,9/2, {i.e.^ the energies at the back of the 

Zel'dovich spiral). The result is 

g/gc-l~ l-E/E,' ^ ' 

where cl and cr are constants unimportant to our overall discussion. Choosing two points 
equally spaced away from gc, call these gi = gc — A^', g2 = gc + ^9^ and their corresponding 
energies Ei = Ec + AE, E2 = Ec — AE we may use two versions of the approximation in 



Eq. (19) to write 

gc ^osc CrEc 

The important point to take away from this analysis is that the width of the rearrangement 
region, i.e. the range in g over which the rearrangement occurs, is ^^ ^ b/losc- An analogous 
result to Eq. (20) is briefiy discussed in Ref. [4J in the context of exotic atoms. There, the 
width is stated to be ^ 6/a^, where as is the Bohr radius. We see that our Eq. (20) is 
consistent with the result for exotic atoms, in that the width of the rearrangement region is 
of the order of the range of the potential over the characteristic length of the problem. 

In Fig. 7, we numerically verify our analytical expression, viz., Eq. (20), by plotting the 
lowest two branches of the FSW for decreasing values of the range, 6, of the potential. It 
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is evident that as 6 ^ 0, the level rearrangements curves evolve to a series of staircase 
functions, which is entirely expected given the collective results of Equations (18) and (20). 
This staircase property of the spectrum has also been discussed in Refs. [4lfT7] in the context 
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FIG. 7: (Color online) Level rearrangement spectrum for several different values of the range, b 
for the FSW. The solid line (orange online) is b/iosc — 0.01, The dashed line (black online) is 
b/f-osc — 0.001, the dot-dot-dash line (red online) b/iosc — 0.0001 and the dot-dash line (green 
online) is b/iosc — 0.00001. Inset: A magnification of the data in the main figure near ^c,o further 
illustrating how the level rearrangements evolve to the staircase profile as 6 ^ 

of the quantum defect of atomic physics, but using an entirely different approach to the one 
presented here. Note that in the inset to Fig. 7, all four curves intersect at a common point, 
namely, dX g = Qc which corresponds to the back (z.e.^ |(/)| = tt) of the Zel'dovich spiral. 

There are two noteworty points to be taken from this staircase like behaviour. The first 
is that any other panel of the spectrum, e.g. gi^ g2 etc. in Fig. 5 (provided b/£osc <^ l), 
can be obtained by simply applying the scaling transformation, Eq. (10), to the data in the 
go region . In addition, the staircase property of the level rearrangements as 6 ^ is not 
specific to the FSW, which implies that for any short-range two-body potential, the E vs. g 
curves will exhibit the same scaling symmetry provided the critical values, gc^n are properly 
scaled as 6 ^ 0. It is also important to realize that even the staircase level rearrangements 
are mapped onto the universal E vs. a spectrum, just as with the other potentials listed in 
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Eq. (7) with h/(.osc <C 1 in the right panel of Figure 1. 

IV. CONCLUSIONS 

In this paper, we have examined the two-body problem of ultra-cold harmonically trapped 
interacting atoms and its relation to the Zel'dovich effect. We have shown, through our con- 
struction of the "Zel'dovich spiral" , that the universal spectrum in terms of the scattering 
length is exactly equivalent to the Zel'dovich effect. This non-trivial observation has been 
used to motivate further experimental studies in order to provide additional data for the 
E vs. a spectrum, which may then be used to establish the first direct experimental obser- 
avtion of the Zel'dovich effect. Finally, we have shown that in the 6^0 limit, the level 
rearrangement spectrum exhibits an exact scaling symmetry, which has until now, gone un- 
noticed. The exact mapping means that the entire E vs. g spectrum (and therefore the E 
vs. a spectrum) may be obtained solely from knowledge of the g^ region as 6 ^ 0. 
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